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Page 5, column 1, line 32: Insert the following after equa- 
tion (19): This expression is similar to Serein's form 
of the Bobyleff-Forsythe equation (Serrin, J.: Math- 
ematical Principles of Classical Fluid Mechanics. 
Encyclopedia of Physics, Vol. VIII/1, Fluid Dynam- 
ics I, S. Flugge and C. Truesdell, eds., Springer- 
Verlag, 1959, p. 251): 

1 o _ 

D:D = - 2 -© +diva 

where 

1 (duj dui, 'l , - _ Du 
D= l[^ + ‘5^J and a_ Dt 

For incompressible, viscous flow: 

Page 7, column 1, line 16: Insert the following after equa- 
tion (34): Equation 34 is a special case of a more 
general result obtained by Rayleigh, who extended 
Helmholz’s work (Lamb, H.: Hydrodynamics. Sixth 
ed. Dover, 1945, p. 617). From a perturbation analy- 
sis, Rayleigh obtained the minimum dissipation 
condition: 


Page 8, column 2, line 39: The phrase “the x -derivatives” 
should be ’‘the ^-derivatives”. 

Page 9, column 1, line 4: The denominator of the second 
term in the second line of equation (48) should be 
squared, not cubed; the calculation used the correct 
expression. (See also other corrections to equa- 
tion (48) below.) 

Pages 9-10. The symbols x and y should be replaced with 
the symbols £ and q in several equations (note that u 
and v are velocity components in the airfoil coordi- 
nate frame). Correct as follows: 


d(u + iv) _( d(u + iv)V dzV^i^ 

— ( dz JWUJ 

r -./■ \ 2 
. 2a 2 e MX 2ik$a sin(a) z 2 

|_( z + $) 3 (z + 8) 2 (z 2 -l) 



V 2 V = VH , H is any scalar 

Rayleigh stated that this condition was exclusive; 
flows that did not satisfy the condition were not char- 
acterized by minimum dissipation. Onsager later 
extended the minimum condition to all steady, linear 
processes that satisfied the second law of thermody- 
namics. Rayleigh used examples of parallel channel 
flow (parabolic profile) and rotational flow, as dis- 
cussed by Lamb, to illustrate his minimum condition. 



du _ dv 3v _ 3u 

5n“d^’ 3q 8 ? ( 50 ) 




Plage 10, column 1, lines 5-14: Insert a period after "zero” 
in line 5, and delete "and results in a net drag . . . 

(51) which is drag times freestream velocity.” 


| <MA V = (DU + W blc )/l = -ypU 3 cC DE 


Ay 



Cde = 



Page 10, column 2, line 26: Insert the following: For all 
the examples of incompressible and fully inotational 
flow, die net drag is zero and the dissipation is due 
entirely to the work of the moving-surface boundary- 
layer control system (W B lc )• This quantity repre- 
sents the minimum power expended by an ideal 
airfoil The drag of an actual airfoil can be related to 
this quantity by defining a boundary layer control 
efficiency: 


T lBLC = 


( DU+W ^U„ 


equal lift, velocity 


Page 9, column 2, line 9: After the word “points.” insert 
the following: An inverse cubic approximation for 
dissipation was used in the leading edge (LE) and 
trailing edge (TE) regions, with a grid spacing of 0.02 
in the LE region (-2.19 < ^ < -1.81, -19 < tj < 0.19) 
and in the TE region (1.81 <£ <2.19,-19 <n < 
0.19). The dissipation approximation was integrated 
analytically in the £ direction and numerically in the 
T1 direction. 


Page 10, column 2, line 32: Insert the following: The 
effective drag coefficient Cde as presented in the text 
represents a power loss coefficient, indicating 
boundary later control power only for the airfoil in 
fully inotational flow and drag plus boundary layer 
control power for airfoils in rotational flow. 
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NOMENCLATURE 


a B 

surface of body subjected to drag force 

L 

airfoil length 

Ay 

surface bounding control volume 

L 

wall length 

a 

circle radius for streamlined-airfoil flow 

1 

channel-wall running length 

a 0 

speed of sound in undisturbed fluid 

M 

sum of terms in mass equation 

B 

constant of integration 

Mj,M k 

Mach numbers in j and k directions 

C 

constant of integration 

n 

unit vector normal to surface 

CdE 

effective-drag coefficient 

P T 

torque power absorbed by rotating cylinder 

c 

airfoil chord 

P 

pressure 

c 

wave speed 

P 

time average (mean) of pressure 

D 

airfoil drag 

P* 

fluctuation in time of pressure 

djk 

deformation-rate tensor 

q 

heat-flow vector 

ee 

unit vector in the 0 direction, cylindrical 
coordinates 

R 

airfoil Reynolds number 

F 

sum of terms in momentum equation 

R c 

critical Reynolds number 

f.f'.f" 

similarity function and its derivatives 

Rs 

jet Reynolds number based on jet half-width 

g.g’.g" 

auxiliary similarity function and its derivatives 

r 

radial coordinate, cylindrical coordinate 
system 

h 

jet nozzle height 

i*LE 

airfoil leading-edge radius 

Lliv 

cross-stream integral of velocity divergence 

n^E 

airfoil trailing-edge radius 

It 

cross-stream integral of total dissipation 

ro 

cylinder radius (vortex-core radius) 

ii 

cross-stream integral of mean pressure 

S 

system entropy 

h 

cross-stream integral of squared mean vorticity 

s 

specific entropy (entropy per unit mass) 

13 

cross-stream integral of the mean square of 
fluctuating vorticity 

T 

temperature 

K 

bulk viscosity coefficient 

t 

time 

ks 

vortex strength factor 

t/c 

airfoil thickness-to-chord ratio 

kx 

coefficient of thermal conductivity 

U 

freestream velocity vector 

L 

torque on a cylinder 

U m 

maximum jet velocity 


iii 
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u,v,w 

velocity components, Cartesian coordinate 

e 

angular coordinate, cylindrical coordinate 


system 


system 

Ui,Uj,U k 

tensor velocity components 


second coefficient of viscosity 

U0 

tangentail velocity component 

P- 

first coefficient of viscosity 

u 

velocity vector 

P 

density 

V 

control volume 

0 

fluid stress 

V c 

channel flow rate 

Ore 

shear stress, cylindrical vortex 

W 

complex potential 

x'jk 

viscous stress tensor 

Wblc 

boundary -layer-control system power 

<D 

viscous dissipation function 

w 

channel width 

O* 

dimensionless dissipation function 

x,y,z 

position components, Cartesian coordinate 


modified dissipation function 


system 


modified dissipation function for incompress- 

Xi,Xj,X k 

tensor position components 


ible flow 

z 

complex position, not transformed 


radiative component of modified dissipation 



function 

a 

Lagrange multiplier, channel flow 

♦ 

unsteady, complex disturbance stream function 

a 

angle of attack 

9 

complex amplitude of the disturbance stream 

a 

wave number 


function 

r 

circulation, cylindrical vortex 

9i 

imaginary component of the amplitude of the 



disturbance stream function 

AS 

rate of change in the system entropy 

9r 

real component of the amplitude of the distur- 

A e S 

contributions to AS from external sources 


bance stream function 

AjS 

contributions to AS from internal sources 

Q 

cylinder rotation rate 

5 

circle offset distance 

to 

complex frequency 

Sj k 

Kronecker delta 

CO 

vorticity 

§1/2 

jet half-width 

CO 

time average (mean) of vorticity 

£ 

disturbance amplitude 

CO' 

fluctuation in time of vorticity 


complex position, transformed 

COj 

imaginary component of complex frequency 

n 

similarity parameter 

Wjk 

rate-of-rotation tensor 


iv 


SUMMARY 


The objective of this report is to examine recent developments in the field of nonequilibrium 
thermodynamics associated with viscous flows, and to relate these developments to the understanding 
of specific phenomena in aerodynamics and aeroacoustics. A key element of the nonequilibrium theory 
is the principle of minimum entropy production rate for steady dissipative processes near equilibrium, 
and variational calculus is used to apply this principle to several examples of viscous flow. The paper 
begins with a review of nonequilibrium thermodynamics and its role in fluid motion. Several 
formulations are presented of the local entropy production rate and the local energy dissipation rate, 
two quantities that are of central importance to the theory. These expressions and the principle of 
minimum entropy production rate for steady viscous flows are used to identify parallel-wall channel 
flow and irrotational flow as having minimally dissipative velocity distributions. Features of 
irrotational, steady, viscous flow near an airfoil, such as the effect of trailing-edge radius on 
circulation, are also found to be compatible with the minimum principle. Finally, the minimum 
principle is used to interpret the stability of an initially laminar, parallel shear flow with respect to 
infinitesimal and finite amplitude disturbances; the results are consistent with experiment and with 
linearized hydrodynamic stability theory. These results suggest that a thermodynamic approach may be 
useful in unifying the understanding of many diverse phenomena in aerodynamics and aeroacoustics. 


INTRODUCTION 

The crucial role thermodynamics plays in fluid 
mechanics and related disciplines such as aerodynamics is 
indicated by the importance of the perfect gas law and the 
energy equation, which, together with the principles of 
mass and momentum conservation, form the basis of the 
theory of fluid motion. The gas law and the energy equa- 
tion can be formulated for a system of fluid particles in 
thermodynamic equilibrium, a system devoid of spatial or 
temporal variations in material properties. For compress- 
ible flows, such as flows containing shock waves, the 
equilibrium approximation is still accurate locally because 
of the large number of particles in the smallest regions of 
interest. 

Spatial and temporal variations in the fluid properties 
may become significant to the thermodynamic description 
of the system when the variations in thermal energy are 
independent of the variations in fluid pressure, such as in 
low-speed flow, or when the variations in fluid properties 
with time are large. When the flow is steady, the theory of 
linear nonequilibrium thermodynamics leads to several 
useful concepts, including the principle of minimum 
entropy production for near-equilibrium processes and the 
subsequent formulation of variational methods consistent 
with this principle. In certain unsteady flows that are 
thermodynamically far from equilibrium, such as Benard 
convection, the heat transfer and entropy production can 
be maximized. 

Variational principles associated with energy dissipa- 
tion of a viscous fluid have been applied with varying 


degrees of success. In the field of hydrodynamics, 
Helmholz proved that flow with negligible inertia forces 
is characterized by a minimum of viscous dissipation 
(ref. 1). Tadmor defined a minimum entropy principle in 
the context of stable numerical solutions of the gas 
dynamics equations (ref. 2). For thermal convection and 
turbulent flow in channels, Malkus proposed a theory 
based on the arguments that the smallest turbulent scales 
are marginally stable, and the mean and turbulent charac- 
teristics of the flow are determined by the condition that 
the energy dissipation is maximized (refs. 3,4). Reynolds 
and Tiederman contended that extrema of the dissipation 
function are not generally compatible with the conserva- 
tion laws, and they found poor agreement between an 
eddy viscosity simulation of turbulent channel flow and 
Malkus’ prediction of the association of maximum dissi- 
pation rate with neutral stability of the mean flow (ref. 5). 
Gage et al. argued against the existence of a general 
thermokinetic variational principle (ref. 6). 

These efforts raise questions about (1) whether a gen- 
eral variational principle is applicable to problems in fluid 
mechanics, (2) which fluid property is to be extremized, 
and (3) whether a maximum or minimum of the property 
should be expected. These issues can be clarified by the 
review that follows of the results of studies of nonequilib- 
rium thermodynamics. For a full discussion of nonequilib- 
rium thermodynamics, the reader is referred to the 
treatises by Yourgrau et al. (ref. 7) and de Groot and 
Mazur (ref. 8). 

Thermodynamic equilibrium, characterized by the 
absence of gradients of state variables or kinematic 



variables, corresponds to a state of zero entropy produc- 
tion and maximum entropy. Linear nonequilibrium 
processes, studied by Onsager (ref. 9), are associated with 
nonzero gradients of system properties and are character- 
ized by minimum entropy production or energy dissipa- 
tion. This characterization specifies that the entropy 
production can be evaluated at every point in the fluid, 
and that the spatial distribution of this quantity is such that 
the volume integral of the entropy production is smaller 
for the observed spatial distribution than for any other 
possible distribution. The volume integral of the entropy 
production is also described as stationary since a small 
deviation from the observed distribution results in a vari- 
ation of the volume integral that is on the order of the 
square of the deviation. The linear description applies if 
the generalized thermodynamic fluxes (such as heat flow 
or fluid deformation) are linearly related to the general- 
ized thermodynamic forces (such as temperature gradient 
or viscous stress). Many near-equilibrium processes are 
accurately characterized as linear. Minimum entropy 
production for these processes can be associated with the 
stability of stationary or time-invariant states in the pres- 
ence of fluctuations in one or more of the thermodynamic 
forces or fluxes. 

Prigogine noted that the entropy production charac- 
teristics for nonequilibrium processes that are nonlinear or 
far from equilibrium can be quite different than those 
observed for processes close to equilibrium (ref. 10). Non- 
linearity is generally associated with the generation of 
chaotic states from initially determinate conditions, such 
as the transition from laminar flow to turbulence. 
Prigogine (ref. 11) and Prigogine and Stengers (ref. 12) 
also demonstrated that nonlinearity and processes that are 
far from equilibrium may give rise to global organization 
from an initially random field, and they characterized 
such processes as the evolution of dissipative structures. 
For example, the problem of convective heat transfer can 
be characterized by the maximum heat transfer and 
entropy production that is associated, in some circum- 
stances, with the appearance of orderly Benard convection 
cells (refs. 13-15). Lugt commented on the evolution of 
discrete vortices in shear flows as another example of 
these dissipative structures (ref. 16). Ziegler proposed that 
rates of entropy production and energy dissipation are 
maximized for processes far from equilibrium (ref. 17). 
These results suggest that maximum energy dissipation 
may be associated with some processes far from equilib- 
rium, or with the transition of a system between equilib- 
rium states. 

The dissipation function in its conventional form 
consists of terms formed from velocity-field gradients. 
Vortices were cited by Lugt in reference 16 as examples 


of dissipative structures in fluid flow, and pressure-wave 
source regions may be examples as well. The first goal of 
this study is to formulate the dissipation function in terms 
directly related to these phenomena. In this context, we 
will derive in the next section an alternative expression of 
the dissipation function; its terms are associated with the 
processes of irreversible compression, vorticity, and con- 
vected wave propagation of density disturbances. 

Many references cited in this discussion refer to the 
steady motion of a viscous fluid as an example of a linear, 
nonequilibrium process for which entropy production 
should be minimized. The second goal of this study is to 
present some specific examples of this concept. Finally, 
we will compare the hydrodynamic and thermodynamic 
stability characteristics of a viscous, parallel shear flow. 

DISSIPATION FUNCTIONS FOR A NEWTONIAN 
FLUID 

Entropy Change in a Fluid Particle System 

Yourgrau et al. studied the rate of change of entropy 
in an open system that allows the transfer of mass and 
energy across system boundaries (ref. 7). For such sys- 
tems, the rate of change of the system entropy, AS, is 
composed of contributions from external and internal 
sources: 

AS = A e S + AjS (1) 

The first term, A e S, represents the entropy exchange 
across the system boundaries by mass or heat transfer. It 
may be positive or negative depending on whether 
entropy is transferred in or out of the system. The second 
term, A|S, represents the entropy production within the 
system boundaries caused by irreversible processes, such 
as heat conduction and viscous dissipation of energy. 
According to the second law of thermodynamics, it is 
always positive. 

If we consider only viscous and thermal conduction 
effects, we can express equation (1) in terms of integrals 
over the system volume, V, bounded by the system 
surface, A: 


2 




-JJf afc-JJWsi* 


A„S 


jjjih 


VT 
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The dissipation function can be written in Cartesian 
coordinates as 



(7a) 


The viscous dissipation function, O, is the inner product 
of the viscous stress tensor, Tj k , and the deformation-rate 
tensor, dj k . 

The surface integrals may be changed to volume 
integrals by using the divergence theorem. Since the 
equality must hold for arbitrary volumes of integration, 
the integrand of the volume integral satisfies the equality 
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Applying the equation of mass conservation gives 
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(7b) 


Ds q <X> q VT 

p D7 ~ ‘ T T “ 


(3b) 


We now restrict our attention to a Newtonian fluid, in 
which the viscous stress is linearly related to the deforma- 
tion, or strain, rate, 

^ jk = 2pdj k + ^jk == 2pdjj c + K-^— — 8j^ (4) 


where 


K = ?i + — u 
3^ 


(5a) 




auj_ + au jL 

3x k 3*j 


(5b) 


Equation (3b) for a Newtonian fluid can be written as 


Ds O k T V 2 T 
P Dt ” T + x 2 


( 8 ) 


Since shear, dilatation, and heat transfer can occur inde- 
pendently, the second law of thermodynamics requires 
that the coefficients \ i, K, and k^ each be positive. The 
second law of thermodynamics can be expressed either 
globally (an isolated system must exhibit zero or positive 
change in entropy) or locally (the viscous dissipation 
function must be positive everywhere and heat flows only 
from hot to cold regions). In the remaining discussion, 
entropy changes caused by heat conduction will be 
neglected, so the entropy production will be considered to 
be related only to the energy dissipation caused by 
viscous effects. 

Alternative Forms of the Dissipation Function 


d jk = 


3u j , du k 

3xk 3xj 


1 ^Uj s /C X 

-33^* (5c) 


and the heat flow is linearly related to the temperature 
gradient, 


q = — kj V T 


(6) 


The dissipation function can also be expressed in 
terms of physical mechanisms that dissipate energy. Using 
Cartesian tensor notation with the summation convention 
on repeated indices, Ziegler (ref. 17) wrote the dissipation 
function of equation (7a) in the form 


® = Xd i jd ick + 2|id JK dj k (9) 
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where d:^ is the fluid-deformation-rate, or strain-rate, 
tensor, defined by equation (5b). Jeffreys (ref. 18) rewrote 
equation (9) as 


. i 3u j 3u k 

O = A, -5 — — 3— - + L1 

oxj dx k 


1 duj du k 

2 a) j k ®jk +2 a Xk a Xj 


( 10 ) 


where 0)j k is the rate-of-rotation tensor, defined by 


0}j k = 


au k 9u j 

3xj dx^ 


(ID 


The last term in equation (10) appears in several 
aeroacoustic wave equations, such as that derived by 
Phillips (ref. 19), Horne and Karamcheti (ref. 20) have 
derived a different formulation of the dissipation function 
by adjoining functions of the conservation equations to 
the dissipation function of equation (10): 

<I* = 0-2pjvF-^|t+MVuJ (12) 

The symbols M and F refer to the sums of terms in the 
conservation equations of mass and momentum, respec- 
tively. They are defined by 
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and 
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(16b) 
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(16c) 


_ dp du k dp du k 9uj 0 U . 
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T i k_P [dx k + dxj J + 5 J kA 'dx j 


(15) 


For a fluid free of sources of mass or momentum, 
M = 0 and F = 0 . 


When terms are canceled as indicated, equation (12) 
becomes 
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(17) 


Expanding the terms of equation (12) individually, 
we obtain 


This expression, which is exact for a Newtonian 
fluid, accounts for three important mechanisms of energy 
dissipation in real flows: (1) irreversible expansion or 
compression of the fluid, (2) the generation and radiation 
of sound or shock waves, and (3) the generation of vortic- 
ity. The modified dissipation function of equation (17) is 
equal to the original dissipation function if the 


4 


conservation equations of mass and momentum are each 
satisfied, so extrema of the volume integral of 
equation (17) should also be compatible with the 
conservation equations. Although equation (17) was 
derived in a manner similar to the method of imposing 
equality constraints on variational problems, the 
conservation equations were not adjoined independently, 
so further consideration of this issue is necessary. 


The radiative terms in equation (17) can be further 
manipulated to yield a convective wave operator acting on 
the density field: 
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►( 18 ) 


These expressions for the dissipation function 
(eqs. (7a), (19)) will be used in the following sections to 
examine several simple flows, to determine if minimally 
dissipative distributions correspond to observed flows. 


INCOMPRESSIBLE-FLOW DISTRIBUTIONS 

The dissipation functions derived in the previous 
section can be used to determine whether simple, steady, 
viscous flows may be characterized by velocity distribu- 
tions that generate stationary or minimum values of the 
net dissipation, as predicted by the theory of linear, 
nonequilibrium thermodynamics. Examples of such flows 
to be considered in the following sections are flow in a 
parallel-wall channel, unconstrained viscous flow, flow 
near a rotating cylinder, and flow near a streamlined 
airfoil. 

Parallel-Wall Channel Flow 


where ao is the speed of sound in an undisturbed fluid, 
and Mj and are the Mach numbers in the j and 
k directions. The first term occurs in Lighthilfs acoustic 
stress tensor (ref. 21), and the second term, a convective 
wave operator for the density, is identical to the convec- 
tive wave operator for acoustical disturbances in shear 
flow (ref. 22). Flows that generate and radiate acoustic 
energy have nonvanishing levels of viscous dissipation in 
the aeroacoustic source region of the flow. (Dissipation in 
the propagation region caused by viscosity and relaxation 
was studied by Meixner, as discussed by de Groot and 
Mazur (ref. 8).) The dissipation function should be posi- 
tive everywhere, although the terms associated with 
outward acoustic radiation are negative, since energy is 
radiated away from, rather than deposited within, the 
aeroacoustic source region. 

For incompressible flow, the dissipation function of 
equation (17) reduces to 

4), =-^V 2 p + |iC0 2 (19) 


We will examine parallel channel flow, or Poiseuille 
flow, as a simple illustration of the method to be used 
throughout this discussion. For this flow, the streamwise 
velocity, u, is a function of the cross-stream coordinate, y, 
only. The viscous dissipation function given by equa- 
tion (7a) reduces to 





2 


( 20 ) 


The net dissipation per unit length per unit width is given 
by 



where 1 and w are the channel-wall running length and 
channel width, respectively. The flow rate per unit width 
through the channel is given by 


The squared-vorticity term, pco 2 , which also appears in 
equation (10), is similar to the quantity enstrophy, which 
is one half of the vorticity squared and is of significance 
in atmospheric flows. Some studies have suggested that 
this quantity is minimized for certain atmospheric vortices 
(ref. 23). The squared-vorticity term is also similar to the 
energy dissipation in homogeneous turbulence, given by 
Hinze (ref. 24) and Tennekes and Lumley (ref. 25) as the 
viscosity times the mean square of the fluctuating 
vorticity. 



( 22 ) 


Possible extrema of the integral in equation (21), sub- 
ject to constant channel flow rate, are found as solutions 
of the Euler-Lagrange equation, 


3f 

3u 



(23) 
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where 


The net dissipation may be expressed in integral form 



f = n 


3u 

By 


+ au 


(24a) 


(24b) 


Here, a is a Lagrange multiplier adjoining the integrand 
of the flow-rate constraint to the extremizing function. For 
this example, evaluation of the Euler-Lagrange equation 
yields the expression 


a = 2\i 



(25) 


The net dissipation is stationary for the parabolic profile 


as 
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Odv 
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= JJJ°( u X’ u y’ u Z’ v X’ v y’ v Z’ w x- w y- w z) dx dy dz 


(29) 


Velocity distributions that generate a stationary or mini- 
mum value of the net dissipation must satisfy the Euler- 
Lagrange equations (ref. 26). For three dependent vari- 
ables (u, v, w) and three independent variables (x, y, z) 
these equations are 


30 d ( d<&\ 
3u 3x(^3u X J 
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(30a) 


u(y) = ^-y 2 +By+c 


(26) 


where B and C are constants of integration. This profile 
corresponds to a minimum of the dissipation function, 
according to Legendre's test (ref. 26), since 


3 2 f 
3u'3u' 


= +2 


(27) 


Unconstrained Viscous Flow 


The preceding example incorporated a constraint on 
the net flow rate. In addition, the equation of mass was 
satisfied by the specified restrictions on the velocity field. 
We now consider general incompressible viscous flows 
that satisfy no constraints other than the integrability and 
differentiability requirements of a variational analysis. 

The dissipation function of equation (7a), when sim- 
plified for incompressible flow, becomes 



where the subscripts x, y, and z refer to partial differentia- 
tion. We will consider various circumstances under which 
the integral of the dissipation function (the net dissipation) 
over a flow region may be stationary. 
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No boundary conditions or other constraints have been 
imposed, although such constraints may be adjoined with 
Lagrange multipliers, as in parallel channel flow. The 
terms in equation (30a) are evaluated as 



(31a) 


30 3u d f 30 ^ 3 2 u 

3u x ^3x 3x(3u x ) ^ 3x 2 


(31b) 


30 

3Uy 


( dv 3u ^ 

3 

30 1 

= 2p 

( 3 2 v + dV 

v 3x 3y J 

; 9y 

3u y J 

^3x3y 3y 2 j 


(31c) 


30 ( 3u 3w \ 

3u z ^\3z 3 x j 


3 ( 30 
3z^3u z 




+ 


3 2 w ^ 

3z3 x j 
(3 Id) 


Summing these terms and rearranging give the Euler- 
Lagrange equation for u, 
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2nV 2 u + 2Hj|-(V-u) = 0 (32) 

Adding the Euler-Lagrange equations for the three com- 
ponents gives the vector equation 

2p.V 2 u + 2|i.V(V • u) = 0 (33) 

The presence of the second term in equation (33) indicates 
that velocity distributions that generate stationary values 
of the net dissipation do not necessarily satisfy the con- 
servation equations, a concern raised by Gage et al. 
(ref. 6) and Reynolds and Tiederman (ref. 5). Extrema 
could be constrained to satisfy the conservation equations 
by adjoining the equations to the dissipation function with 
influence functions. Instead, we impose conservation of 
mass by setting the second term equal to zero. The net 
dissipation is then stationary for velocity distributions that 
satisfy 

V 2 u = 0 (34) 

This condition is equivalent to the vanishing of the 
viscous-force vector term in the momentum equation for 
incompressible flow, 

P^ = - v P + H y2 a (35) 

This condition is also equivalent to a zero value for the 
curl of the vorticity field throughout the flow; this can be 
seen with the vector identity 

V 2 u = V( V -u)-Vx(VxD) (36) 

The dissipation function is zero for the case of uni- 
form rotation, since the fluid deformation is zero. The net 
dissipation is also stationary for velocity distributions that 
are irrotational. The viscous force on a volume element of 
fluid is proportional to the divergence of the stress and is 
zero for irrotational flow. However, the viscous stress, 
and hence the dissipation, are not generally zero for 
regions in a viscous fluid where the flow is irrotational. 
An analogous situation is found in the dissipative process 
of steady heat flow within a conducting solid, where the 
irrotational heat flux, q , is proportional to the gradient of 
the scalar temperature, T. In the following sections, we 
examine two familiar flows with regard to their descrip- 
tions as irreversible thermodynamic processes. 


Flow Near a Rotating Cylinder 

The incompressible flow near a hollow rotating 
cylinder is a special case of flow with regions of solid 
rotational and irrotational motion. It is described in cylin- 
drical coordinates by 

u = u e (r)e e (37) 

For this axisymmetric tangential flow, the viscous dissipa- 
tion function of equation (7a) reduces to 

o, r)=11 ^_ia) 2 (38, 

and the Euler-Lagrange equation for the net dissipation 
reduces to the Euler equation, 

r 2dfuj, +r due o <39, 

dr 2 dr 

The two solutions are 

ue«r,r _1 (40) 

These two solutions represent respectively the rotational 
core inside the cylinder and the potential region of a 
columnar vortex outside the cylinder. Both solutions rep- 
resent minimum dissipation conditions, as in the channel 
flow. 

Integrated contributions to the overall dissipation rate 
in the field of a cylinder with radius ro and circulation T 
are given in table 1 . As can be seen, the net dissipation in 
the rotational core is zero, while the irrotational outer flow 
has a net positive dissipation. The tangential acceleration 


Table 1. Velocity distribution and net dissipation in a 
steady vortex 


Region 

u e( r ) 

JJordrdG 

0 < r < r 0 

Tr 

0 

+ o 

2ktq 



r 

ur 2 

r 0 < r < oo 

+ 2nr 

+ 2 
1*0 
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is zero throughout the flow, except at the boundary at tq. 
A steady tangential stress is required at this location 
(provided by a thin rotating cylinder, for example) to 
maintain a steady flow. The net dissipation in the outer 
region is identical to the torque power absorbed by the 
rotating cylinder, since the shear stress at the cylinder wall 
is given by 


ars(ro) = H 


9 uq 

9r 


= _iiL 
r ' r 0 nT 0 


(41) 


of air around airfoils and streamlined bodies. A 
few decades ago the potential theory was con- 
sidered a purely mathematical discipline. Later 
it was understood that it might be a valuable 
guide for understanding of the general laws 
which govern the lift, moment, induced drag, 
etc. of airfoils. However, as aerodynamical 
design became more refined, it was found that 
the quantitative predictions of the theory are of 
a remarkable accuracy in spite of its idealized 
conception. 


The cylinder torque is then given by 
2n 

L = -|o^(r 0 )r^ d6 = 2nr 
0 


This, combined with the cylinder rotation rate, 

Q=5M.i T 

r 0 2uTq 


The high degree of correspondence between associ- 
ated viscous and inviscid flows implies that for many 
flows, viscous effects are confined to a thin layer adjacent 
(42) to a moving surface, and that large regions of the real flow 
sufficiently far from boundary surfaces are irrotational. In 
this section we evaluate the distribution of the viscous 
dissipation function in the vicinity of a streamlined airfoil. 

The airfoil flows investigated in this study were gen- 
erated by mapping an offset circle in the z-plane to a 
' ' streamlined airfoil shape in the i^-plane with the confor- 
mal transformation 


gives the torque power absorbed by the rotating cylinder, 

M r 2 

= (44) 

71 r 0 

In this example, the stationary dissipation is associ- 
ated with either uniform or irrotational motion, as in 
unconstrained viscous flow. Yates and Donaldson 
(ref. 27) and Greene (ref. 28) have also considered energy 
dissipation in vortices in their analyses of induced drag. 


; = z + - (45) 

Z 

The radius of the circle in the z-plane is a, centered at 
z = -8. A sharp trailing edge is generated if a = 1 + 8. A 
round trailing edge is generated if a > 1 + 8. The circle in 
figure 1 maps to an airfoil with a chord length of approx- 
imately 4. The complex potential, W, caused by a doublet, 
a variable-strength vortex, and a uniform stream of speed 
U inclined at angle a can be derived from the expression 
given by Abbot and von Doenhoff (ref. 30), 


Flow Near a Streamlined Airfoil 

Another important example of viscous irrotational 
flow is the low-speed flow of a real fluid over a stream- 
lined object such as an airfoil. At a sufficiently high 
Reynolds number, in the absence of separation, the 
viscous-flow effects are confined to a thin boundary layer 
near the surface, a fact first utilized by Prandtl in his 
boundary layer theory. As von Karman noted in 1941 
(ref. 29): 

One of the most surprising developments in 
modern fluid mechanics is the successful 
application of the theory of potential motion of 
ideal incompressible fluids to the actual flow 


(z + 8)’ 


W = ui (z + 8)e -iot + — - e ia + 2ik s a sin(a) 


x In 


(z + S)e 


-IOC 


(46) 


The vortex strength, kg, is equal to 1 when the separation 
streamline departs from the airfoil trailing edge, and 
kg = -1 for stagnation at the airfoil leading edge. For 
sharp trailing-edge airfoils in an inviscid fluid, the Kutta 
condition requires that kg = 1 to avoid infinite velocities 
at the trailing edge. The velocity and the x-derivatives of 
the velocity components are given by 
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Figure 1. (a) Cylindrical flow in z coordinate system . (b) Airfoil in transformed (Q coordinate system. 


dW dz 

u + iv = 

dz dq 


3u _ 3v 3v _ 3u 

3y 3x ’ 3y 3x 


(50) 


_ia a 2 e ia 2ik s asin(a) 

z 2 ) 

(z + 8) 2 (z + 8) 

U 2 -iJ 


(47) 


3(u + iv) 
3x 


/ d(u + iv) 
dz 


dz V 3^ 
d<; Jv3x J 


so the dissipation function becomes 


<D(x,y) = 4|i 



(51) 


= U 


' 2a 2 e ia N 
(z + 5) 3 


2iksasin(a) 
(z + 8) 3 


' z 2 ^ 


U 


a 2 <- ia 

— i(Y d C 

e 7T + 


p-') 

2iksasin(a) 


(z + 8) 2 ( Z + S) 


2z J 


2z 3 




(48) 


For an incompressible, two-dimensional flow, the dissipa- 
tion function reduces to 



+ 


f 3v 3u ^ 
^dx dy, 


2 


(49) 


If the flow is irrotational as well as incompressible and 
two-dimensional, then 


The dissipation can be nondimensionalized using the 
chord, the viscosity, and the speed: 


O* = 



(52) 


The dissipation function thus obtained was evaluated 
on a square grid with a spacing of 0.10 between grid 
points. Figure 2 illustrates the streamlines and contours of 
constant dissipation in the vicinity of a Joukouski-type 
airfoil, with boundary layers excluded. The irrotational 
flow model corresponds to a high degree to the associated 
real flow, provided the Reynolds number is sufficiently 
large and no regions of separated flow are present. The 
irrotational model is even more accurate if the real airfoil 
is provided with a moving-surface boundary-layer-control 
system that can adjust the tangential velocity distribution 
over the airfoil surface to match the distribution predicted 
by the inviscid theory. This method was used by Yates to 
study unsteady airfoil flow (ref. 31). The thin boundary 
layer itself is a good approximation of such a control 
system. 
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where 



& 
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Figure 2. Symmetrical airfoil in transformed coordinate system, 
t/c = 0.1, a = 0.1 rad; streamlines are dashed lines and contours 
of dimensionless dissipation are solid lines. 


As mentioned in the section on unconstrained viscous 
flow, the viscous force per unit volume in the fluid, which 
is proportional to the divergence of the stress, is zero 
throughout the irrotational region. The viscous stress, 
however, is not zero and results in a net drag on the airfoil 
even if the boundary layer is absent (assuming the no-slip 
condition is satisfied with a moving-surface boundary- 
layer control). This drag can be computed by integrating 
the streamwise component of the dot product of the 
viscous stress with the outward normal vector around a 
control surface such as the airfoil surface. The drag can 
also be found by using the energy equation; the volume 
integral of energy dissipation can be set equal to the drag 
power, which is drag times freestream velocity. 
Lagerstrom derived an expression for the drag power DU 
(ref. 32), 


u-U l2 


DU = JJp j-Lu nda- JJ G -(u-Oj-nda 

A v A v 

+ JJJ a : ^dv + JJ ° u ■ nda (53) 


A v is the surface bounding the control volume 
Ab is the surface of the body subjected to drag force 

U is the freestream velocity 

V is the control volume 

o is the fluid stress 

The terms on the right-hand side of equation (53) 
account for, in order, (1) net kinetic energy transfer out of 
the control volume, from wakes or trailing vortices, 
(2) stress work on the control surfaces, (3) stress work 
within the control volume, and (4) stress work on the sur- 
face of the body subjected to the drag force. The last term 
is zero if the body surface is motionless, but nonzero if the 
surface is moving. For the case of two-dimensional flow 
with negligible wakes or vortex, this expression can be 
simplified to 


DU + Wblc 



(pV • u + <J>)dv 


V 


(54) 


where Wgoc is equal to the negative of the fourth term 
of equation (53). The right-hand side of equation (54) is 
equal to the third term of equation (53). For incompress- 
ible flow, the effective drag power per unit length, 
(dU + Wblc)/l> is the surface integral of the 
dissipation. 


J <MA V = (DU + W BLC )/L = jpU 3 c C DE 
A v 



(55) 


where c is the airfoil chord and L is the airfoil length. 
This gives an expression for the effective drag coefficient. 



where the Reynolds number is R = pUc/p. The variation 
of Coe R f° r two zero-incidence airfoils that have 
different thickness ratios is compared in figure 3 with that 
of a flat plate that has boundary layers on both sides. 
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(a) 


r TE /r LE 



R 

Figure 3. Zero-incidence effective-drag coefficient vs. Reynolds 
number for a flat-plate boundary layer ( two sides) and for 
symmetrical airfoils (boundary-layer drag excluded on airfoils). 



a, rad 


Figure 4. Variation of effective-drag coefficient with angle of 
attack, R = 1000. 

The variation of the effective-drag coefficient with 
angle of attack is shown in figure 4 for airfoil thickness- 
to-chord ratios of 0.1 and 0.2. The effective-drag coeffi- 
cient is minimized at zero incidence and increases as the 
square of the incidence angle, with the higher rate of 
increase associated with the airfoil that has the lower 
thickness-to-chord ratio. The region of highest flow curva- 
ture and dissipation is situated at the airfoil leading edge, 
as seen in figure 2. 

Both the trailing-edge velocity and the net dissipation 
obtain infinite values if the airfoil circulation is not suffi- 
cient to maintain the separation streamline at the airfoil 
trailing edge. The association between net dissipation and 
trailing-edge streamline location can be examined further 



(e) ( - 0.05 

Figure 5 . Airfoil section profiles with rounded trailing edges; 
rjE/c is (a) 0.0068, (b) 0.0045, (c) 0.0027, (d) 0.0013, (e) 0.0007 . 

by considering the hypothetical case of continuously 
adjustable circulation around airfoils with rounded trailing 
edges (and with moving-surface boundary-layer treat- 
ment, as already discussed in this section). 

Five 16.4%-thick, rounded-trailing-edge airfoil pro- 
files were considered for this study, with trailing-edge 
radius-to-chord ratios of 0.0068, 0.0045, 0.0027, 0.0013, 
and 0.0007 (fig. 5). The geometrical parameters for these 
profiles are given in table 2. 

The net dissipation at an incidence of 0.1 rad was 
found by evaluating the velocity derivatives analytically 
and integrating numerically, as described earlier in this 
section. The variation in effective drag coefficient with 
circulation for these airfoils at R = 1000 is shown in fig- 
ure 6. For the elliptical airfoil (case 1), the net dissipation 
is minimized for zero circulation. The circulation corre- 
sponding to minimum net dissipation increases as the 
trailing-edge radius, rpE» is reduced relative to the 
leading-edge radius, i*le (fig- 6). 
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Table 2. Parameters for rounded-trailing-edge airfoil 


Parameter 

Case 

1 

2 

3 

4 

5 

radius, a 

1.18 

1.178 

1.172 

1.165 

1.160 

offset, 8 

0.000 

0.033 

0.0613 

0.090 

0.105 

LE a radius 

0.0275 

0.0368 

0.0441 

0.0518 

0.0551 

TE b radius 

0.0275 

0.0184 

0.0110 

0.0052 

0.0029 

r TE /r LE 

1.0 

0.500 

0.249 

0.101 

0.0521 

r TE /c 

0.0068 

0.0045 

0.0027 

0.0013 

0.0007 


a LE = leading edge 
b TE = trailing edge 



Figure 6. Drag coefficient vs. circulation coefficient; R = 1000, 
t/c = 0 . 164, a = 0.1 rad. 


Streamlines and contours of constant dimensionless 
dissipation for full circulation (trailing-edge separation) 
and for circulation corresponding to minimum dissipation 
for the third case are shown in figure 7. For minimum 
dissipation conditions, the rear separation point moves 
rearward and the separation-stream angle decreases, as the 
trailing-edge radius decreases. Figure 8 illustrates how the 
kg needed for minimum dissipation varies with ryg/r^g; 
this variation is nearly linear on a semilog plot. These 
results suggest that the familiar circulation condition for 
airfoils with sharp trailing edges can be approached in the 
limit as the trailing-edge radius vanishes. The results also 
suggest that the minimum dissipation condition associated 
with thermodynamic stability is sufficient to uniquely 
determine the circulation for airfoils with finite trailing- 
edge radii. 

Attempts to compare these results with experiment 
must take into account the presence of the boundary layer 
and the zones of separation near the rounded trailing edge. 
For an airfoil with boundary layers, terms associated with 
rTE/ r LE should have an effect of order 1/R, while terms 
associated with the ratio of trailing-edge radius to 
boundary-layer thickness will have a stronger effect, of 
order 1/VR. Thwaites notes that the Kutta condition, 
which is required to avoid a singularity at the trailing edge 
of an airfoil in an inviscid fluid and which uniquely 
determines the value of circulation, has no counterpart in 
the case of viscous flow. Thwaites suggests a circulation 
condition for airfoils in viscous flow as the condition that 
the rear separation streamlines originate from the trailing 
edge at an angle that lies between the extensions of the 
upper and lower trailing-edge surfaces (ref. 33). 

These results do not demonstrate that the entire field 
(irrotational regions and boundary layers) of a real viscous 
flow satisfies a condition of stationary dissipation; how- 
ever, the value of the dissipation function for any irrota- 
tional region will be much less than the value in rotational 
regions. Consequently, a condition of minimum net dissi- 
pation in a real flow may correspond to the establishment 
of irrotational flow to the largest extent possible. 

This analysis can be extended to include the effects of 
boundary layers, if desired. The association between 
thermodynamic stability and minimum net dissipation can 
also be extended to examine problems of hydrodynamic 
stability, as will be seen in the next section. 
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Figure 7. Streamlines and dissipation contours in transformed coordinate system, full circulation (upper) and minimum drag circulation 
(lower); rjg/c = 0.0027, k s for minimum dissipation = 0.570. 








r TE /r LE 

Figure 8. Circulation coefficient for minimum drag vs. /TEALS' 
a = 0.1, t/c= 0.164. 


Viscous laminar flow may also be considered an open 
thermodynamic system. Energy transfers within the flow 
in the forms of heat conduction and energy dissipation 
through viscous stress are governed by the first and sec- 
ond laws of thermodynamics, as discussed in the intro- 
duction. If gradients in mean velocity exist within the 
flow, the system is not in thermodynamic equilibrium and 
the entropy production is nonnegligible. Under certain 
restrictions, the second law of thermodynamics implies 
that the thermodynamic state is stable when the entropy 
production is minimized (ref. 7). This thermodynamic 
stability condition has been used to examine the hydro- 
dynamic stability of fluid motions such as pipe Poiseuille 
flow (ref. 37), thermal (Benard) convection (ref. 38), and 
transonic flow (ref. 39). However, the understanding of 
the extent of the correspondence between hydrodynamic 
and thermodynamic stability of fluid motion is not com- 
plete and would benefit from further clarification. To this 
end, a simple numerical illustration was generated from 
the results published by Chun and Schwartz of a stability 
analysis of a laminar wall jet (ref. 35). This exercise is 
described in the next section. 


HYDRODYNAMIC AND THERMODYNAMIC 
STABILITY CHARACTERISTICS OF A LAMINAR 
WALL JET 

The initial stages of the transition of laminar flow to 
turbulence have been extensively studied by means of lin- 
earized hydrodynamic stability theory. This analysis has 
been performed for many two-dimensional, parallel shear 
flows, such as the boundary layer, channel flow, the free 
jet (ref. 34), and the wall jet (refs. 35,36). The analysis has 
also been extended to nonparallel flows, three- 
dimensional flows, and compressible flows. In general, 
the theory predicts that many laminar flows are stable 
below a critical Reynolds number, R c , whose value 
depends strongly on the velocity profile. Below R c , dis- 
turbances are damped. Above R c , disturbances are ampli- 
fied exponentially in either space or time at a rate that 
depends on the disturbance frequency. At any R > R c , 
there is a frequency corresponding to the maximum 
amplification rate; this frequency is comparable to the 
experimentally observed frequency of unforced distur- 
bances in the shear flow. Experimental disturbances in 
shear flows with initially low turbulence levels exhibit 
good agreement with linearized hydrodynamic stability 
theory in regions where the disturbances are amplified 
through several orders of magnitude. At high disturbance 
levels, the flow may undergo secondary instabilities and a 
transition to turbulence. 


Wall-Jet Hydrodynamic Stability 

In their stability analysis of the laminar wall jet, Chun 
and Schwarz predicted a critical Reynolds number of 58, 
with amplified disturbances in the form of a double vortex 
row similar to a Karman vortex street adjacent to the wall. 
They considered only temporally amplified disturbances, 
with real wave numbers and complex frequencies. Tsuji 
et al. confirmed the Chun and Schwartz results with a 
study that also examined spatially amplified disturbances 
(real frequencies and complex wave numbers) and veri- 
fied the prediction with experiment (ref. 36). 

The mean velocity field considered in these studies, 
as well as in the present investigation, is a parallel flow 
derived from the Glauert similarity solution (ref. 40): 

u = 3.1749U m f (57a) 

v = 0 (57b) 

y/5 1/2 =0.3955ti (57c) 


1 + g + g‘ ) n _i J30 

— - ^ — +V3tan 1 g (58a) 

0-g) L(2+g)J 
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f’=2gg', g'=I(l-g 3 j 


(58b) 


4 


f = f(Tl), g = g(rj) (58c) 

Profiles of the velocity and its first derivative are plotted 
in figure 9. The disturbance stream function takes the 
form of a convecting wave, 

< H x >y> t) = <p(y)exp[ia(x - ct)] (59) 



Figure 10. Stability characteristics of the ideal wall jet (ret 36). 



Figure 9. Profiles of the wall-jet mean velocity and its first 
derivative. 

The disturbance stream function satisfies the Orr- 
Somerfeld equation, the boundary conditions at the wall, 
and the conditions far from the wall. The boundary condi- 
tions, which define an eigenvalue problem, allow the 
complex frequency, co = ac, to be determined as a func- 
tion of the Reynolds number, R§, and the real wave 
number, a. The plot by Tsuji and Morikawa (ref. 36) of 
the imaginary component 0 )j of the complex frequency is 
reproduced in figure 10. Chun and Schwarz presented the 
real and imaginary components of the disturbance stream 
function for three conditions, 

1. a = 1, Rg = 30 (stable) 

2. a = 1, Rg = 58.2 (neutral) 


described by Chun and Schwarz. A typical disturbance 
streamline pattern for the wall jet is shown in figure 12, 
which shows the double-row vortex structure characteris- 
tic of the wall jet. 


Wall- Jet Thermodynamic Stability 

The disturbance velocity field was evaluated for these 
three conditions in order to compute the viscous dissipa- 
tion function, which for incompressible, two-dimensional 
flow reduces to 


o=n 



+ 2 


f dv'l 

dy 


+ 


< av + au > 

dy j 


2 


(60) 


The mean-velocity derivatives, evaluated from the Glauert 
solution, are 


3u 

dy 

= 3. 1749U m f" 

(61a) 

du 

3v 

3v 

(61b) 

= — = — = O 

3x 

dy 

3x 


where 

f”=2(g') 2 + 2gg'\ g"=-g 2 g' (62) 

The velocity derivatives of the disturbance, evaluated 
from the disturbance stream function at t = 0, are 


3. a = 1 , R§ = 150 (unstable) 

These conditions are indicated in figure 10 and the stream 
function components are presented in figure 11. The 
stream function components were extended from y/Si /2 = 
2.4 to y/Sj/2 = 4.0 with a decaying exponential, as 


3u' 

3x 


3v 

dy 


e[cx(p r cos(ax) + acpi sin(ax)] 


(63a) 


3i f 

dy 



cos(ax) - 


d 2( Pi 

3y 2 


sin(ax) 


(63b) 
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Figure 11. Real and imaginary components of the disturbance stream function; a=1,R s = 30, 58.2, 150 (ref. 35). 



ax 


Figure 12 . Disturbance streamlines for an ideal wall jet; R = 200, 
c = 2 (ref. 36). 


— = e[a 2 (p r cos(ax) - a 2 <pj sin(ax) (63c) 

3x L 


This expression was then integrated numerically over the 
region 0 < x/5i/ 2 < 4 and 0 < y /5 l/2 < 2n t0 evaluate the 
net dissipation. For this purpose, the integral was divided 
into three terms. 
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(65) 


Here, e is a small parameter that is proportional to the 
amplitude of the disturbance. Negative values of 8 are not 
permitted because they would violate the boundary condi- 
tions on the disturbance stream function. Errors associated 
with inaccuracies in the linearized equations for small but 
finite amplitude disturbances are neglected here. 

The velocity derivatives of the mean and disturbed 
flow were used to evaluate the total viscous dissipation. 
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(64) 


For a given condition of Reynolds number and wave 
number, the value of the net dissipation should vary 
quadratically with the disturbance amplitude and have 
positive curvature everywhere. The variation of net dissi- 
pation for the three Reynolds number conditions 
examined here is shown in figure 13. The variations of 
dissipation for negative values of disturbance amplitude 
are shown as dotted lines, since these variations violate 
the boundary conditions. If negative amplitudes were 
permitted, the corresponding vortical disturbances would 
rotate in the sense opposite to the actual disturbances. 
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Figure 13 shows that for R = 30, the condition asso- 
ciated with stable laminar flow, the value of the net dissi- 
pation increases with increasing disturbance amplitude. 
For R = 30, the minimum value of net dissipation is gen- 
erated for vanishing disturbance amplitude, suggesting 
that this flow is stable thermodynamically as well as 
hydrodynamically. For R = 58, the case of neutral hydro- 
dynamic stability, the curve of net dissipation has a hori- 
zontal tangent at the origin, suggesting that the flow is 
thermodynamically neutral for infinitesimal disturbances 
but stable for large disturbances. For R = 150, the net dis- 
sipation decreases as the disturbance amplitude increases 
from 0, until a minimum value of the net dissipation is 
reached at e = 0.005. This suggests that the hydrodynami- 
cally unstable mean flow is also thermodynamically 
unstable, although the disturbed flow is thermodynami- 
cally stable for the value of e = 0.005, at which the net 
dissipation is minimized. This situation may correspond to 
the nonlinear limit cycle that is sometimes observed in the 
real flow. Figure 14 shows a flow visualization of an ini- 
tially laminar wall jet at R = 600; the large disturbances 
form a stable vortex array for many cycles in the down- 
stream direction. Similar patterns form in the downstream 
wakes of cylinders or other bluff bodies (ref. 41). Such 
vortex arrays in an inviscid fluid show good agreement 
between predictions of stable motion and experimental 
results. Vortex arrays in a viscous fluid can also be con- 
sidered from the standpoint of thermodynamic stability. 
This approach may contribute to the understanding of the 
discrete tone stabilization of separation-layer and 


boundary-layer control. In the next section, measurements 
of the velocity field of an unstable, initially laminar wall 
jet are used to estimate the components of the dissipation 
function in various regions of an unstable jet. 

Dissipation in an Unsteady Wall Jet 

To further consider dissipation effects in a real flow, 
it is convenient to separate the dissipation function into 
mean and fluctuating components. The components of 
pressure and vorticity are 

p = p + p\ co = co + go' (66) 


so the dissipation function of equation (22) becomes 

<X>I =-^iv 2 (p + p') + ^i(co^-co , ) 2 (67) 

and the mean dissipation function is 

<j»I =-y-V 2 p + ^(co) 2 +(co') 2 (68) 


Figure 14 illustrates the periodic vortex motion of a two- 
dimensional wall-jet flow. This jet flow developed from 
an initially parabolic profile rather than a self-similar pro- 
file, which was considered in the preceding section. 



Figure 13 . Net dissipation integrated over 0 < x/S-j/g < 2k and 0 < y/$i/2 < 4.0 vs. disturbance amplitude; R = 30, 58.2, 150. 
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where 





=ii +i 2 +13 


(69) 



Figure 14. Wail-jet phase-averaged velocity field and flow 
visualization; jet height = 0.508 cm, jet width = 10.16 cm, wall 
length = 3.81 cm, maximum velocity of the parabolic exit profile = 
13.85 m/s, tone frequency = 600 Hz. a) Phase-averaged flow 
visualization, b) velocity vectors relative to convecting vortices, 
c) vorticity contours. 

Woolley and Karamcheti showed that the stability charac- 
teristics of nonparallel shear flows are closely related to 
those of parallel shear flows (ref. 42). The short wall was 
associated with the generation of tones by the jet; this tone 
generation significantly reduced the small-disturbance 
region by effectively forcing the initial jet region. How- 
ever, the observed tone frequencies were approximately 
the same as those of the jet without forcing. The condi- 
tions for this jet are U m = 13.85 m/s and L/h = 7.5, where 
U m is the maximum velocity of the parabolic exit profile, 
L is the wall length, and h is the nozzle height. Fig- 
ure 14(a) shows a phase-averaged schlieren visualization 
record obtained by lightly heating the subsonic nozzle 
flow. Figure 14(b) shows the phase-averaged velocity 
field, relative to the convecting vortices, and figure 14(c) 
depicts the corresponding vorticity field obtained with a 
central-differenced curl of the velocity field. The vorticity 
is normalized with respect to the maximum exit velocity 
and the nozzle height. We measured the velocity with a 
single x-wire velocity probe by sampling the probe output 
at regular phase intervals, which were determined by a 
pressure transducer fixed in the wall (ref. 43). The mea- 
surements were further processed to obtain estimates of 
the mean and fluctuating components of the dissipation 
field. 

Relative contributions to the dissipative structure of 
the wall jet were estimated by integrating the field vari- 
ables in the cross-stream direction. We define 


The Laplacian of the pressure was computed from the 
divergence of the incompressible, viscous, two- 
dimensional momentum equation, 


"V 2 p 
p 




+ 4 


du A 

dy 


d v 
dx 


+ 2 


3u 

dy> 


(71) 


To estimate the measurement error, the mean square of 
the divergence of the velocity was computed: 


^div - 


h i 2 r r 3 u + 3v^ 
U m J i{dx^ d y> 


2 


d 


(72) 


This integral should vanish if the flow is two-dimensional 
and incompressible, and if measurements of u and v are 
accurate. 

Figure 15 shows the distributions of I], I2, I3, Idiv» 
and Ij for the jet conditions described. The measurement 
error indicated by Idiv is negligible only upstream of the 
wall trailing edge, where the pressure Laplacian term, Ij , 
is also negligible. The square of the mean vorticity, I2 , is 
nearly constant in the small-disturbance region, 0 < x/h < 
2.5, beyond which it gradually decreases. Within the 
large-disturbance region, 2.5 < x/h < 7.5, the net dissipa- 
tion gradually increases with streamwise distance. The 
mean square of fluctuating vorticity, I3, is negligible in 
the initial region, then steadily increases in the large- 
disturbance wall region. The decrease in Ij beginning at 
x/h = 2.5 is consistent with the decrease in net dissipation 
associated with hydrodynamically unstable flow that was 
noted in the preceding section. 
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Figure 15. Distribution of integrated dissipation components. 

These measurements demonstrate a method of 
directly measuring dissipative terms in an unsteady flow. 
This method can be extended to experimentally determine 
a relationship between overall dissipation and variable 
parameters such as forcing frequency, and to experimen- 
tally search for dissipation extrema. 


CONCLUDING REMARKS 

Many features of steady viscous fluid motion, such as 
the large extent of irrotational flow near streamlined 
bodies, are consistent with a thermodynamic interpretation 
based on the minimum dissipation of energy or production 
of entropy. Observations of this agreement have been 
made by various researchers over the last century; 
however, recent developments in the field of nonequilib- 
rium thermodynamics justify further investigation along 
these lines. In this report, an exact equation was derived 
for the dissipation function of a homogeneous, isotropic, 
Newtonian fluid; the equation’s terms are associated with 
irreversible compression or expansion, wave radiation, 
and the square of the vorticity. Simple flows, such as the 
incompressible channel flow and the cylindrical vortex, 
were identified as minimally dissipative distributions. A 
correspondence between hydrodynamic and thermo- 
dynamic stability is suggested by calculations and mea- 
surements of an initially laminar wall jet. These results 
suggest that the theory of nonequilibrium thermodynamics 
can serve to unify the understanding of many diverse phe- 
nomena in aerodynamics and aeroacoustics. 

National Aeronautics and Space Administration 
Ames Research Center 
Moffett Field, California 94035-1000 
June 27, 1990 
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